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THE GEOMETRY OF GRAVITATIONAL TENSING 
MAGNIFICATION 

AMIR BABAK AAZAMI AND MARCUS C. WERNER 


Abstract. We present a definition of unsigned magnification in gravi¬ 
tational lensing valid on arbitrary convex normal neighborhoods of time 
oriented Lorentzian manifolds. This definition is a function defined at 
any two points along a null geodesic that lie in a convex normal neigh¬ 
borhood, and foregoes the usual notions of lens and source planes in 
gravitational lensing. Rather, it makes essential use of the van Vleck 
determinant, which we present via the exponential map, and Ethering- 
ton’s definition of luminosity distance for arbitrary spacetimes. We then 
specialize our definition to spacetimes, like Schwarzschild’s, in which the 
lens is compact and isolated, and show that our magnification function is 
monotonically increasing along any geodesic contained within a convex 
normal neighborhood. 


1. Introduction 

Flux is an important observable in gravitational lensing. Suppose F is the 
observed flux of a light source, which has been increased due to the grav¬ 
itational focusing by an intervening massive lensing object, and that the 
hypothetical flux in the absence of the lens would be Fq. Then the magni¬ 
fication factor due to gravitational lensing is given by 

F 

F = (1) 

ignoring image orientation. This magnification factor is not regarded as 
an observable because the intrinsic flux Fq without the lens is unknown in 
general. The recent report of the hrst strongly lensed type la supernova [10] 
proves an exception to this rule because such supernovae are standardizable 
candles. Moreover, it is clear that many more examples of strongly lensed 
type la supernovae will be found by upcoming surveys, so cases of observable 
lensing magnification will become more common. 

From a theoretical point of view, the magnification factor ^ is usually con¬ 
sidered in terms of the quasi-Newtonian approximation for gravitational 
lensing (see, e.g., [12, 7]), which treats light rays as piecewise straight lines 
in Euclidean space. However, the proper arena for gravitational lensing is, of 
course, the Lorentzian spacetime geometry of General Relativity. Thus, it is 
desirable to generalize the definition of magnification to a spacetime setting 
and better understand its geometrical meaning. The aim of this article is 
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to do this by reinterpreting the classic definition of luminosity distance in 
spacetimes due to Etherington [2] in terms of lensing magnification and the 
van Vleck determinant [14], which we express here in terms of the exponen¬ 
tial map. Indeed, we are not the first to promote the use of the exponential 
map in gravitational lensing; see, e.g., the viewpoint offered in [11], For the 
spacetime view of gravitational lensing in general, see, e.g., [ 6 ]. 

In our understanding of the van Vleck determinant, we have been greatly 
aided by the two very comprehensive treatments [15, 9]. Indeed, most of 
Sections 3 and 4 below, on Synge’s world function and the van Vleck de¬ 
terminant, can be found in [15] and [9], with one exception, however, that 
warrants their inclusion here; namely, our use of the exponential map to 
compute the van Vleck determinant. For this reason, our notation, and 
several of our proofs, are noticeably different from most of the existing liter¬ 
ature. Our definition of magnification (Definition 3) and subsequent focusing 
theorem (Theorem 1) are to be found in Section 5. 

2. Overview of the exponential map and normal coordinates 

Let (M, g) be a connected time oriented four-dimensional Lorentzian mani¬ 
fold, with g having signature (—, -|-, +). Let be a convex normal neigh¬ 

borhood of p € M, that is, a neighborhood any two points q,q' of which are 
connected by a unique distance-minimizing geodesic aqq/ lying entirely in 
though there may well be other geodesics between q and q' that leave 
and then reenter the aqq/ are usually referred to as “radial geodesics,” 
and we will adopt this terminology henceforth. Furthermore, one can ar¬ 
range it so that for each (7 G is contained in the normal neighborhood 

of q provided by the exponential map at g; that such a ^ exists at every 
point on a Lorentzian manifold is proved, e.g., in [4, pp. 133-136]. Thus 
normal coordinates (x*) centered at any point p £ ^ cover all of ‘if, and 
we will in fact describe Synge’s “world function” below primarily in terms 
of such coordinates. Next, “aqq” will always denote the unique radial ge¬ 
odesic in ^ from q to q', ‘^Oi'^qi{t)” will denote its tangent vector at the 
point Oqqi (t) , and we henceforth adopt the convention of parametrizing our 
radial geodesics to run for unit affine parameter: aqqi: [ 0 , 1 ] — ‘if, with 
Oiqq'{0) = q and aqq'{l) = q'. Also, we adopt the convention of writing both 
a curve, and its coordinate representation in a coordinate chart, using the 
same symbol. Finally, the Einstein summation convention will be assumed 
throughout, with indices labeled 0 , 1 , 2 ,3. 

Because of the essential role played for us by the exponential map, we now 
briefly review some of its properties. Thus, fix p £ and recall that any 
point q £ is given hy q = yviX) some unique vector V £ 

TpM, where yv is the unique geodesic starting at p in the direction V, and 
where exp^ denotes the exponential map at p. Since exp^ is a diffeomorphism 
from a neighborhood of 0 G TpM to any choice of orthonormal basis 
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{Eq\p, ..., E'sip} C TpM provides us with “normal” coordinates (x®) defined 
with respect to that basis. Indeed, expressing any X S TpM as X = X''Ei\p, 
the diffeomorphism E: TpM —M” sending X i—^ composes 

with exp“^ to give 

q = ^v{l) ^ = V^Ei\p A {V\...,V^) ;= x{q) . (2) 

^ V 

(ai*) coordinates of q 


Thus x'^{q) = V\ and because geodesics jv have the scaling property 
7v(i) = 7tvW whenever either side is defined, the geodesic jy in the coor¬ 
dinates (x®) is given by 

jv(t) = (3) 


Of course, yy: [0,1] —> ^ must be our radial geodesic Opg: [0,1] —> 
Furthermore, the normal coordinate basis vectors {d/dx^,... ,d/dx^} de¬ 
fined with respect to the (x®) satisfy, by construction, 


A 

5x® 


= Ei 


P) 


hence gij{p) = diag (— 1 , 1 , 1 , 1 ) (i.e., they are “normal” at 
(3), it also follows that r®^;,(p) = 0, hence also di\p{gjk) = 
[5, Prop. 33, p. 73]. 


p). Because of 
0 ; consult, e.g.. 


Now we use normal coordinates at p to define “quasi-normal” coordinates 
at any other point g € as follows. Let {FloIp; ■ ■ ■ ^E^lp} C TpM denote 
the orthonormal basis with respect to which the normal coordinates (x®) 
at p are defined. Given any other point q G let Jj denote the unique 
Jacobi field along the radial geodesic Upg: [0,1] —> ^ satisfying Ji(0) = 0 
and J{(0) = Ei\p^ where “J'” denotes the covariant derivative of Jj along 
apq. Observe that {Jo(l), ... , J 3 (l)} C TgM will be linearly independent 
provided that q is not conjugate to p along Upg. Indeed this must be so, 
since ^ is contained in the normal coordinate chart centered at p, no point 
in ^ is conjugate to p along a radial geodesic through p (this is an important 
point to which we will retnrn later; consult, e.g., [5, Prop. 10, p. 271]). Thus 
{Jo(l),..., ^ 3 ( 1 )} C TgM is a basis, though it need not be orthonormal. In 
any case, use it to define “quasi-normal” coordinates (x®) centered at g € 
via the exponential map exp^ at q, in the same manner as in ( 2 ) above. 
Then, by construction. 


A 

5x® 


Q 


= Ji{^) 


for each i = 0,..., 3. In fact, each 


A 

clx® 


A 

5x® 


Q 


g 


( 4 ) 
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as follows. Set q = expp(y) and let (t>Ei\p\v £ TyiTpM) denote the vector 
canonically identified with V € TpM (explicitly, C°°{TpM) —)• M 

is defined by / e-)- Then a standard resnlt relating 

exponential maps and Jacobi fields (see, e.g., [5, Prop. 1, p. 217]) says that 

dviexpp){(j)Ei\p\v) = Ti(l), (5) 

where dv'(expp): Tv(TpM) — )• TqM is the differential of exp^ at V. Finally, 
observe that if g' = expq(lF) and W = VF*Ji(l) G TqM, then the radial 
geodesic : [0,1] —)■ ^ is given in the coordinates (x*) by aqqi{t) = 
{tW^,... ,tW^). We will make nse of these coordinates in Section 4 below. 

3. The world function and its gradient 

Though, in fact, we are already in a position to define the van Vleck deter¬ 
minant, in this section we first define Synge’s world function [13] and state 
some of its properties. We do this primarily to connect our work with previ¬ 
ous presentations of the van Vleck determinant appearing in the literature, 
in particular [9], from which this section derives. 

Definition 1. With notation as above, define the world function 

1 

, cr{p,q) := - g{a'pq,apq) dt. (6) 

Note that when q = p, the radial geodesic app is just the constant curve 
Oipp{t) = p, so that a'ppit) = 0 G TpM and cr{p,p) = 0. In fact, cr is a 
smooth function (see [5, Lemma 9, p. 131]) that is invariant under a change 
of parametrization, in the following sense: if we instead parametrize our 
radial geodesic apq to have domain [a, b] and define 

, . b — a . / X , 

J 9Wpq,apq)dt, 

then this function is invariant under any geodesic reparametrization of 
Upq: [a, b] —)• ‘tf. Note that if L{apq) denotes the length of apq: [0,1] —> 
then 

{ T{apqfi/2 apq spacelike, 

0 apq null, 

—L{apqfif2 apq timelike. 

Moving on, observe that (6) is very similar to the energy function of any 
curve segment a: [0,1] —> 

1 

F;a:[0,l] —, Ea{t) ■■= - g{a'{u),a'{u))du 

^ Jo 

(see [5, pp. 288ff.]). Indeed, let x : [0,1] x {—5, J) —^ be a smooth variation 
of apq through geodesics starting at p; i.e., each curve x(-,v): [0,1] —> ‘if 
is a geodesic starting at p. Denoting the tangent vectors of each x{-,v) 
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by dx{d/dt\(^. ;= Xt{-,v), we define the energy Ex - {—S,6) —> M of the 
variation x of Opq to be the smooth function 


Ex{v) 


1 

2 



g{xt{-,v),xt{-,v))dt. 


Then it is straightforward to verify that a{p, ■) o a;(l, •): (—5, 5) —>■ M sat¬ 
isfies 


aip, ■) ox{l,-) = Ex- (7) 

If our parametrization domain were [a, b], then (7) would be a{p, •)oa;(l, •) = 
{b — a)Ex- This fact plays an important role in the proof of the following 
fundamental result regarding the world function a{p, ■): ^ —)■ M with initial 
point p fixed. 


Lemma 1. Fix p G and let apq'. [0, 1] —^ he the radial geodesic in ^ 
from p to an arbitrary point g G In normal coordinates (x*) centered at 
p, the gradient grad(j(p, •) := K at q is 


Kq — apg(l) 

Moreover, K satisfies 


xfiq) 


A 




( 8 ) 


g{K,K) = 2a{p,-). 


(9) 


Proof. Throughout this proof, we follow the notation and terminology in [5, 
pp. 215ff.]. If x{q) = {x^{q),... ,x^{q)) in normal coordinates (x*) centered 
at p, then 

= {tx^{Q),---,tx^{q)). ( 10 ) 

Now let Wq = di\q G TqM be arbitrary and consider the variation 
X : [0,1] X (—(5, (5) — > ^ defined in the normal coordinates (x*) by 

x{t,v) := (t{x^{q) + vW^),... ,t{x^{q) + vW^)). 

This is well-defined and contained in ‘if for 5 small enough, with all lon¬ 
gitudinal curves x{-,v): [0,1] —being geodesics hxed at p, base curve 
a;(-,0): [0,1] —>■ ‘tf our original geodesic apq in (10), and with variation 
field V: [0,1] — > TM given by 



Because the variation is fixed at p, 

y(0) = OG TpM, 

while at the opposite end 


V{1) = Wq€ TqM. 
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Note that F is a Jacobi field because all longitudinal curves are geodesics. 
To derive (8), we now consider the differential of the composition 

a{p,-)ox{l,-): (-6,6) —>R , v ^ a{p,x{l,v)), 

and find that at u = 0, 


do{(^{Pr) o®(l,-)) 


I dv 




= g{Kg,Wg). 

The significance of this can be seen once we realize that, thanks to (7), 

{a{p,-) ox{l,-)){v) = E^{v), 

and the variational properties of the latter are well known. Indeed, for the 
variation x, 

dEx 


dv 


gi^oipq , 17 ^ 


= 9{OLpq{l),Wq), 


(for the full form of £1^(0), see [5, Prop. 39, p. 289]). Noting that our choice 
of Wq ^ TqM was arbitrary, it follows that we must have 

A 

(9a:* 


Kq = a'pqil) = x\q) 


We point out here that if Opq: [a, b] —then the right-hand side of (8) 
would be scaled by a factor (b — a). Working in normal coordinates centered 
at p so that we can rely on (8), we evaluate g{K, K) at q £‘if and obtain 

g{K,K){q) = g{a'pq{q),a'pq{q)) = 2a{p,q), 

where we note that the function t g{a'pq{t),a'pq{t)) is constant along the 
geodesic apq. □ 

The following fact, whose proof we omit, now follows as a consequence. 


Corollary 1. Fix p £ and let K = grad(T(p, •). Then Kq ^ 0 for all 
q £‘f — {p}, while Kp = 0. Furthermore, K satisfies 

VkK = K, 

hence its integral curves admit parametrizations as geodesics. 


Before investigating the integral curves of the vector field K, we ask: In 
Lemma 1, what if we had wanted to fix the endpoint q instead, so that we 
would be working with a{-, q ): —> R? Not surprisingly, the answer is the 
same, as we state here without proof. 


Corollary 2. Fix q £ ‘if and let Opq: [0,1] —>■ ‘if be the radial geodesic in ‘f 
from an arbitrary point p £ ‘if to q. In normal or quasi-normal coordinates 
(y*) centered at q, the gradient giada{-,q) := K at p is 


d 


Kp = -a'pqiO) = fiP)-^ 


( 11 ) 
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Moreover, K satisfies 


g{K,K) 


2a{-,q). 


( 12 ) 


We close this section by describing the integral curves oi K = grad a{p, •). 
Because Kp = 0 by Corollary 1, the integral curve of K through p is simply 
the constant curve at p. Now pick any other point g € ^ and let 7 : I —?■ ^ 
denote the maximal integral curve of K starting at q, where I C M is a 
connected open interval containing 0. In terms of normal coordinates (x*) 
centered at p, let 7 (f) = ( 7 °(t),..., 7 ^(f)) and x{q) = {x^{q),... ,x^{q)). 
Then 


lit) = fi\t) 


A. 


7h) 


= = Yit) 


A. 


7(b 


where the last equality arises via (8) in Lemma 1 . Since 7(0) = 
{x^{q),..., x"^{q)) , it follows that in normal coordinates (x*) centered 
at p, the integral curve of K starting at q is 


lit) = (e*x°(g),...,e‘x^(g)). 


Hence the maximal interval I necessarily contains (—oo,0] and 


lim 7 (f) = ( 0 ,..., 0 ) = p, 

t^ — OO 

though p ^ 7 (/). Of course, S/yj' = 7 ' by Corollary 1 , hence 7 (f) admits a 
reparametrization as a geodesic. It is straightforward to verify that one such 
reparametrization is given by the diffeomorphism h: I —)• h{I) C (0,+oo) 
defined by /i(f) = e* ;= s, in terms of which the reparametrized curve 
7 ( 5 ) := {y oh~^){s) is given, once again in normal coordinates (x*) centered 
at p, by 

7(s) = (sx°(q),...,sx^(q)) , ^(s) = —lllllfl. 


Observe that 7 is certainly defined on (0,1] = h((—oo,0]), on which it must 
coincide with the radial geodesic ( 0 ,1] —>■ and that 7(1) = q. 


4. The van Vleck determinant 

The quantity known as the van Vleck determinant [14] arises via a simple 
application of the Jacobian function of a smooth mapping (see [5, p. 196]). 
Recall that, given smooth oriented manifolds M and N with corresponding 
volume forms ujm and wat, and a smooth map fi: M — N, the Jacobian 
function of (p is the smooth function ^{4>) € C°°{M) such that 

(P*{un) = (13) 

where p* is the pullback of p. Now let p be the diffeomorphism 
exp^: expp^('^) C TpM —)■ The volume element on in terms 

of normal coordinates (x*) centered at p, is det [gy] dx^ A ••• A dx^, 
while the volume element on exp“^(‘^) C TpM, which we identify via the 
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difFeomorphism E in (2) with an open subset in is simply dx^ A • • ■ Adx^. 
Then (13) becomes 

exp* det [gij] dx^ A • • • A dx^'^ = ^(exp^) dx^ A • • • A dx^. 

For any q £ ^ with corresponding coordinates x{q), observe that 
dexpp{d/dx^\x(q)) = d/dx''\q, hence 

^(expp)(g) = -det [gij {q)], (14) 

where, by identifying q = expp(y) € with its unique vector V € TpM, we 
can regard ^(exp^) as a smooth function on rather than on TpM. Note 
that, using the quasi-normal coordinates of Section 2, and (4) in particular, 
the right-hand side of (14) is also given by 


^-det[5((Ji(l), Jj(l))]. 


(15) 


Indeed, recalling that 7 \/(t) = 7t\/(l) whenever either side is defined, it 
is straightforward to show that for all t > 0 at which 'jv is defined, the 
differential dtvi^^Pp)- Ttv{TpM) —^ satisfies 


(t>Ei\p\tv ^ dtv{e^Vp){(t>Ei\p\tv) 


Ji{t) 


whose derivation is similar to that in [5, Prop. 1, p. 217]. Then the Jacobian 
function of exp^ along 'yv satisfies 


^(expp)(7y(i)) 


v'-det [g{Ji{t),Jj{t))] 


Definition 2 (van Vleck determinant in a convex normal neighborhood). 
Let be a convex normal neighborhood of a Lorentzian manifold {M,g). 
The van Vleck determinant A: ^ x ^ —> M is the function defined by 


A{p, q) 


1 

^(expp)(g) 


(16) 


Let us make hve remarks about this dehnition. (1) Recall that, because 
we have restricted ourselves to a convex normal neighborhood any two 
points p,q £ ^ are connected by a unique geodesic segment in if, so there 
is no ambiguity as to the choice of geodesic. (2) We reemphasize that using 
the convex normal neighborhood ^ allows us to write ^(expp), which is 
dehned on TpM, as a smooth function on ^ itself. (3) Bearing (14) in mind, 
note that, although each orthonormal basis for TpM determines a normal 
coordinate system centered at p, A(p, q) can be evaluated with respect to any 
such coordinates: orthonormal bases are related by an orthogonal change 
of basis matrix, hence A(p, q) remains unchanged. (4) Because p = expp(O) 
and doexpp: TQ{TpM) — TpM sends (f)x\o i-A X, the van Vleck determinant 
satishes 


A{p,p) = 1 . 


( 17 ) 
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(5) In Minkowski spacetime which is itself a convex normal neighborhood 
of each of its points, and whose van Vleck determinant we denote henceforth 
by “Aq,” it is easy to verify that 

Ao{p,q) = 1 (18) 

In fact, (16) is not the usual way in which the van Vleck determinant is pre¬ 
sented. Rather, it is defined via the world function a. In [9, 15], for example, 
the van Vleck determinant is defined via the mixed partial derivatives of a: 


dx^dx^ 

Specifically, it is defined to be the following quantity: 

det [-crii{p,q)] 


A{p,q) := - 


V-det [gij{p)] ^-det [gjj{q)] 


(19) 


As is easily verified, this quantity is independent of the coordinates (x*) 
at p and (x*) at q used to compute it. In which case, let us evaluate (19) 
using normal coordinates (x*) centered at p and corresponding quasi-normal 
coordinates (x*) centered at q; doing so, and making use of (4), we will show 
that one obtains simply 

y/-det [gij{q)y ^ 

which is equal to 1 /^ {expp){q) in the same coordinates. To derive (20), 
return to the quantity 

d2. 


a 


dx^dx^ 


{p,q) 


and note that, by Corollary 2, in quasi-normal coordinates (x*) centered at 
q, K = grad( 7 (-, 5 ) satishes 




dx'- 


d 

dx^ 


= x^{p) 


d 


dx^ 


= -«w( 0 )- 


( 21 ) 


j- ^ p 


Now introduce normal coordinates (x*) centered at p, so that we can express 
—apq(O) in two ways: 


_k, . d 


= -«pg(0) = -x'"{q) 


d 


dx^ 


= -x^{q) 


dx^ 


dx^ 


d 


q coordinates 

Together with (21), this yields 
da{-,q) 


p coordinates 


9x® 


= -g-j{p)x [q) 


dx^ 


dx^ 


( 22 ) 
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With q fixed, a{-, q) is a function of p G since both coordinate charts (x*) 
and {xY cover let us “switch charts” and write the left-hand side as 

da{-,q) 


dx^ 

so that (22) can be rewritten as 


dx^ 

da{-,q) 

p dx^ 

p dx^ 


dx'' 


dx^ 


dx^ 


dx^ 

,dx^ 


= -gik{p)x^{q). 


(23) 


9lk{p) 


Now we must take the partial derivative of this with respect to d/dx^, 
evaluated at q. But as we saw in (4), quasi-normal coordinates satisfy 


A 

dx^ 


A 


in which case the application of d/dx^\q to (23) yields simply 


A 

dx^ 

so that, in the end. 


A 

dx^ 


52 


a 


dx^dx^ 


A 


(p-q) 


= -9ik{p), 


= -9lk{p)- 


(24) 


(25) 


With this convenient formulation of the mixed partial derivatives a-j^i of 
a, consider (19) once again. Evaluating it using normal coordinates (x*) 
centered at p and corresponding quasi-normal coordinates (x*) centered at 
q, we obtain, with the help of (25) and (4), simply 

1 

y^-det [gij{q)y 

which is precisely (20). Nevertheless, the reason we have emphasized the 
Jacobian function of the exponential map in Definition 2 is because the 
exponential map gives a precise indication of when q is conjugate to p: 


Lemma 2. For any p £ M and V G TpM, let 'jv denote the geodesic 
starting at p in the direction V. Then the van Vleck determinant A(p, •) is 
unbounded at a point q along 'jv if and only if q is conjugate to p along qy. 

Proof. Let us suppose that q = expp(l/) = 7 y(l). If q is conjugate to p 
along 7 y, then by definition there exists a nontrivial Jacobi field J along 
7 y: [0,1] —)• M satisfying J(0) = J(l) = 0 (see [5, Prop. 10, p. 271]). 
Set J'{0) = W (note that W G TpM must be nonzero, since J is non¬ 
trivial). Then dy(expp)((/)u/|y) = J(l) = 0. Scaling J so that W has 
unit length, it follows from (15) that with respect to any orthonormal ba¬ 
sis for TpM containing W, ^(expp)(g) = 0, hence A(p, •) is unbounded 
at q. Conversely, suppose that A(p, •) is unbounded at q; if the Ji(l)’s 
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were linearly independent, hence a basis for TqM, then the nondegener¬ 
acy of Qq-. TqM X TgM — )• M would ensure that det [ 5 (Ji(l), (1))] ^ 0. 
Since this is not the case, it follows that the Ji(l)’s must be linearly depen¬ 
dent; hence there exists a nontrivial linear combination such that 

X]r=i Jii^) — C* (each ai G M). As any linear combination of Jacobi fields 
is again a Jacobi field, it follows that J := X](Li Ji is a nontrivial Jacobi 
field satisfying J(0) = J(l) = 0 along yy: [0, 1 ] — ^ M. Hence p and q are 
conjugate. □ 

We add immediately that by the Morse Index Theorem for null geodesics 
(see [1, Theorem 10.77, p. 398]), a null geodesic has only finitely many 
eonjugate points. Also, expanding the components of the metric g in normal 
coordinates (x*) centered at p, for points q = 7 y(l) we have 

-det[gijiq)] = 1 - 1 Ric (H, H) + ^(|H|3), 

hence the lowest order expansion of the van Vleck determinant is 

= ^(exp,W(l)) “ ' + 

cf. [15, Eqn. (64)] and [9, Eqn. (7.3)]. 

5. Magnification in gravitational tensing 

Now we come to gravitational leasing in spacetimes. What we need to define 
is a magnification function of two points along a null geodesic in a spacetime, 
the latter modeling a light ray in whose magnification we are interested. As 
this null geodesic is emitted by a light-emitting source, we must therefore 
view our initial point p as lying on the worldline of that source, that is, on 
a smooth, future-pointing timelike curve 7 (r), which, for convenience, we 
take to start at 7 ( 0 ) = p in Let q G ^ and consider the radial geodesic 
apq: [0,1] —> First, we show that the quantity 

da{j{T),q) 

dr r=o 

is nonzero. We will show this via a slight alteration of the proofs of Lemma 
1 and Corollary 2, as follows. In normal coordinates (y*) centered at q, let 
7 *(r) denote the components of 7 (r), so that, in particular, y*(p) = 7 *( 0 ), 
and define a variation x: [0,1] x (—J, J) —> ‘if of Upq by 

x{t,T) := ((1-t)/(r),...,(l-f) 7 ^(T)), (26) 

which is well-defined for 6 small enough. Observe that ® is a variation of 
x{t,0) = Cipq{t) through geodesics, all of which are fixed at the endpoint q. 
The crucial new input here is that the initial transverse curve now coincides 
with the worldline of our light-emitting source: 

x{0,t) = (/(r),...,7^(T)) = 7(r). 
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Now we proceed as we did in Lemma 1 and Corollary 2 and take the deriv¬ 
ative of the composition 

cr(-,g) o ®(0, •): , T a{x{0,T),q) = cr(7(r),g), 

at T = 0. Doing so is made easy once again because 


(cr(-,g) o7)(r) = Ex{t), 


(27) 


with the energy of x satisfying 


dEx 
dr 0 


9iapq,V) 


1 

0 


-fi'(ap5(0),7'(0)) > 0, 


where V{t) denotes the variation field XT-{t,0), which is a Jacobi field sat¬ 
isfying 1/(0) = 7'(0) and 1/(1) = 0, and where g'(Q!pq(0),7'(0)) < 0 because 
Oipq{t) is a future-pointing null geodesic and 7 (t) is a future-pointing timelike 
curve. In other words, the function r i—>■ cr{'y{T),q) necessarily satisfies 


da{'y{T),q) 


dr 


T=0 


-5«(0),f(0)) > 0, 


(28) 


which is nonzero, as desired. Observe that (28) is, in fact, parametrization 
invariant; i.e., if we had worked with a different parametrization of apq, say 
one of the form dtpq: [a,b] — > ^ with apq{a) = p and apq{b) = q, so that 
a'pq{a) = {b — a)~^apg{0), then in fact (27) would be replaced by {cr{-,q) o 
7)('^) — ~ o)Ex{t), so that the right-hand side of (28) would remain 

unchanged. With this established, we are now ready to give a spacetime 
definition of magnification in gravitational lensing. To that end, recall that 
in Euclidean space, an isotropically emitting light source of luminosity L 
has an observed flux E at radial distance r given by 


F 


L 

dvrr^ 


Now, while both L and F are local properties that are meaningful in an 
arbitrary spacetime, the distance r is not. However, one can use L and F to 
construct a distance measure that is meaningful in spacetime, namely, the 
luminosity distance D: 


D{p, q) 


L{p) 

4:'KF{q) 


(29) 


Thus, a naive generalization of the lensing magnification to spacetime follows 
from (1) and (29), 


Kp, q) 


D^{p,q)' 


(30) 


where the subscript 0 denotes “the absence of the lens” — of which we will 
have more to say below. Now, Etherington [2] has shown that the luminosity 
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distance of a light source in an arbitrary spacetime can be written as 


D{p,q) = - 


da{-f{T),q) 


dr 


det [gijjq)] 
det [gij{p)] 


(31) 


provided that this is evaluated in normal coordinates (x*) centered atp, whose 
trajectory in spacetime is given, as above, by a future pointing timelike curve 
7 (t) satisfying 7 ( 0 ) = p. Indeed, in these coordinates |det [( 7 ij(p)]| = 1. 
Therefore, if we combine (14), (16), (30), and (31), we arrive at last at a 
dehnition for the leasing magnihcation in an arbitrary spacetime in terms 
of the van Vleck determinant: 


Definition 3 (Unsigned magnification). Let {M,g) be a spacetime modeling 
a gravitational lens. Let p,q M lie in a convex normal neighborhood of 
M. Then the unsigned magnification of p at g is defined to be 


g{p, q) 



A(p,g) 

^o{p,q)' 


(32) 


where the symbol “0” denotes the spacetime M in the absence of the lens. 


Observe that, because our leasing scenario takes place in a convex normal 
neighborhood of our spacetime, then A(p, q) and {da/dT\o)~^ are both finite 
and positive: finite by (28), and positive by Corollary 2, the latter because 
p will have no conjugate points along any radial geodesic through it. Be¬ 
cause /U > 0 , our definition is that of unsigned magnification; regarding the 
parity of leased images, see, e.g., [12, p. 34], Having said that, note that 
while A{p,q) and da/dT\o are defined on our spacetime M, dcro/drjo and 
Ao(P) 9 ) need not be: “the absence of the lens” may well imply a different 
smooth manifold Mq M, not just a different metric on M. Let us give an 
important example of when this is the case. 

If the mass distribution acting as lens is modeled by, say, the Schwarzschild 
spacetime, then the absence of this lens simply means the limit m —>■ 0 , where 
m is the mass parameter. But the limit m —)■ 0 is Minkowski spacetime M|, 
for which (18) dictates that Aq = 1.^ In fact if we assume that our source 
and observer are “spatially constant,” then we can easily determine the 
quantity dao/drlo as well. Indeed, in terms of global normal coordinates 
(x*) in let us imagine the worldline 7 (r) of our light-emitting source as 
given simply by 

7(r) = (r, 0,0,0). 

In other words, our light source is “spatially constant” in the sense that 
its trajectory is the integral curve of the time orientation d/dx^ starting at 
p = 7 ( 0 ) = 0. Similarly, we imagine our trajectory, as observers, as also 
being spatially constant and hence given by some 

7o(s) = (s,ri,r2,r3). 

^The notion of a “limit of a spacetime” is generally coordinate-dependent; see [3]. 
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This scenario in is now as follows: at r = 0, our light-emitting source 
emits null geodesics in all null directions, precisely one of which will reach us, 
say at 7o(s*) = q. This scenario then repeats itself at every r > 0 thereafter, 
as both source and observer go forward in along their respective integral 
curves of d/dx^. By (28), 


dcro{l{T),q) 


dr 


= -5(aoa(0)’V(0)). 


But since ccog ■ [0,1] —> l^i is given by 


(33) 


aoq{t) = (ts*,tri,tr2,tr3), 


and since it is a null geodesic, it is easily verified that dcJo/drlo is just the 
“spatial separation” between source and observer: 

= + r| + rf := r,. (34) 

0 

Therefore, if we consider the “spatially-scaled” magnification fJ-{p,q)/rg := 
fj,s{p,q), then this would now be a function defined solely on Schwarzschild 
spacetime, and which, bearing (16) and (28) in mind, inspires the following 
definition: 


da'o{'y{T),q) 

dr 


Example (Unsigned magnification; compact, isolated body, with Minkowski 
limit). Let {M,g) be a spacetime modeling a compact, isolated gravitating 
body of mass m whose limit m ^ 0 is Minkowski spacetime. Let 7(t) 
be a future-pointing timelike curve starting at j{0) = p € M and 7y(t) 
any future-pointing null geodesic starting at p in the direction V, such 
that 7v'(l) = q lies in a convex normal neighborhood of p. Then the 
spatially-scaled magnification ps of p at g is 


Ts{p,q) 


1 1 
ff(7'(0),expp^(g))2 ^(expp)(g)’ 


(35) 


where ^(exp^) is the Jacobian function of the exponential map exp^. 


Let us make three remarks regarding this example. (1) Implicit in this defi¬ 
nition is the assumption that the family of metrics g{m) depends smoothly 
on m and has a well defined limit as m —>■ 0; as mentioned above, such an 
assumption is coordinate-dependent in general. (2) Recall our observation 
of the parameter invariance of (28), so that one cannot scale the right-hand 
side of ps by simply rescaling the initial null direction V. (3) In order for 
(32) to be well defined, a Lorentzian manifold Mq representing the situa¬ 
tion of “absence of lens,” as well as points p,q € Mq with which the points 
p,q ^ M can be identified in that absence, must be clearly defined — though 
we do not pursue this here, we point out that one way of making precise the 
identifications of such points is via the method in [3]. Now, as (35) shows, 
when Mq = M|, the determination of the points p, g £ is unnecessary, 
since Aq = 1, and we can scale away the “spatial separation” arising from 
the factor dao/drlQ. Indeed, because we are assuming that the portion of 
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the light ray in which we are interested resides entirely within a (very large) 
convex normal neighborhood, so that there are no conjugate points along 
it, then a further benefit is the following. Before stating this result, let us 
formalize what we mean by a light ray in a spacetime. 

Definition 4 (Light ray). Let {M,g) be a time oriented Lorentzian man¬ 
ifold. A light ray in M is a complete, future-pointing null geodesic 7 sat¬ 
isfying Ric(y,y) > 0 and whose orthogonal complement 7 -*- is everywhere 
integrahle. 

Theorem 1. In a convex normal neighborhood ofp, pLsip,-) is monoton- 
ically increasing along any light ray through p. 

Proof. This follows from the fact that the van Vleck determinant is mono- 
tonically increasing along any light ray in for a proof, see [15]. □ 

We close this paper by pointing out two interesting directions for future work 
that arise within the framework that we have established here. The first is¬ 
sue is to better understand the phenomenon of multiple lensed images in 
this setting. In particular, it would be interesting to see if the occurrence of 
multiple lensed images can occur within a convex normal neighborhood. A 
second related issue is to formulate a spacetime analogue of a magnification 
invariant associated to multiple lensed images. Such invariant sums of the 
signed image magnification are known to occur for certain classes of lens sys¬ 
tems in the standard approximation but there are, as yet, no generalizations 
applicable to spacetimes (see, e.g., [ 8 ] and the references therein). 

Acknowledgements 

This work was supported by the Hakubi Center for Advanced Research, Ky¬ 
oto University, Kyoto, Japan, and the World Premier International Research 
Center Initiative (WPI), MEXT, Japan. The authors thank Eric Poisson 
and Matt Visser for very helpful discussions. 

References 

[ 1 ] J. K. Beem, P. Ehrlich, and K. Easley, Global Lorentzian Geometry, vol. 202, 
CRC Press, 1996. 

[2] I. M. Etherington, Republication of: LX. On the definition of distance in General 
Relativity, General Relativity and Gravitation, 39 (2007), pp. 1055-1067. 

131 R. Geroch, Limits of Spacetimes, Gommunications in Mathematical Physics, 13 
(1969), pp. 180-193. 

[4] N. J. Hicks, Notes on Differential Geometry, Van Nostrand, Princeton, 1965. 

[5] B. O’NEILL, Semi-Riemannian Geometry with Applications to Relativity, vol. 103, 
Academic press, 1983. 

[6] V. Perlick, Ray optics, Fermat’s Principle, and Applications to General Relativity, 
vol. 61, Springer Science & Business Media, 2000. 

[7] A. O. Betters, H. Levine, and J. Wambsganss, Singularity Theory and Gravita¬ 
tional Sensing, vol. 21, Springer Science & Business Media, 2001. 


16 


[8] A. O. Fetters and M. C. Werner, Mathematics of gravitational tensing: multiple 
imaging and magnification, General Relativity and Gravitation, 42 (2010), pp. 2011- 
2046. 

[9] E. Poisson, The motion of point particles in curved spacetime, Living Rev. Relativity, 
7 (2004). 

[10] R. M. Quimby, M. C. Werner, M. Oguri, S. More, A. More, M. Tanaka, 
K. Nomoto, T. J. Moriya, G. Folatelli, K. Maeda, and M. C. Bersten, 
Extraordinary magnification of the ordinary type la supernova PSl-lOafx, The Astro- 
physical Journal Letters, 768 (2013), p. L20. 

[11] P. H. Reimberg and L. R. Abramo, The Jacobi map for gravitational lensing: the 
role of the exponential map, Glassical and Quantum Gravity, 30 (2013), p. 065020. 

[12] P. Schneider, J. Ehlers, and E. E. Falco, Gravitational Lenses, Springer, 1992. 

[13] J. L. Synge, Relativity: the general theory, North-Holland Publishing Company 
Amsterdam, 1960. 

[14] J. H. Van Vleck, The correspondence principle in the statistical interpretation of 
guantum mechanics. Proceedings of the National Academy of Sciences of the United 
States of America, 14 (1928), p. 178. 

[15] M. ViSSER, van Vleck determinants: Geodesic focusing in Lorentzian spacetimes. 
Physical Review D, 47 (1993), p. 2395. 

Kavli IPMU (WPI), UTIAS 

The University of Tokyo 

Kashiwa, Chiba 277-8583, Japan 

E-mail address: aniir.aazaini@ipmu.jp 

Yukawa Institute for theoretical physics 

Kyoto University 

Kyoto 606-8502, Japan 

E-mail address: werner@yukawa.kyoto-u.ac.jp 


